*Replication Instructions for 
*Theorizing Potential Downstream Cultural Consequences of LGBT+ Activism

*Below are instructions for replicating all manuscript and online appendix tables and figures 
*in STATA using the dataset “Social Forces Replication Data.dta”. 

*For Excel charts, use “Social Forces Replication Excel Charts”. 
*Please contact Sam Whitt (swhitt@highpoint.edu) for questions regarding data replication.

*Manuscript Table/Figure Replication

*Figure 2. Changes in LGBT+ Support Over Time Outside Sarajevo

*installing the cibar command
ssc install cibar

cibar gayrights if sarajevocity==0, over1(treatment2) 
graph save g1, replace
cibar supportpride if sarajevocity==0, over1(treatment2) 
graph save g2, replace
cibar closetogay if sarajevocity==0, over1(treatment2) 
graph save g3, replace
cibar dgprogay if sarajevocity==0, over1(treatment2) 
graph save g4, replace
graph combine g1.gph g2.gph g3.gph g4.gph
*Note: Additional Labeling and Formatting Required

 *Table 1. Change in LGBT+ Support Over Time Outside Sarajevo (OLS Regression)
 
reg mmx_progayindex2 i.treatment2 if sarajevocity==0, cluster(municipality)
reg mmx_progayindex2 i.treatment2##c.eth_frac13 voteeu ethnocentric religious  if sarajevocity==0, cluster(municipality)
reg mmx_progayindex2 i.treatment2##c.eth_frac13  voteeu ethnocentric religious hadcovid diedcovid vaccinated Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0, cluster(municipality)
*Note: Additional Labeling and Formatting Required using outreg2 in STATA

*Figure 3. The Impact of Cultural Receptivity on LGBT+ Support (Marginal Effects)

*download the package below to run the <mcp> command
findit gr0056

reg mmx_progayindex2 i.treatment2##c.eth_frac13 if sarajevocity==0, robust
mcp eth_frac13 treatment2, at2(1 2)
graph save g5, replace
reg mmx_progayindex2 c.eth_frac13  voteeu if sarajevocity==0 & yr2021==1, robust
mcp eth_frac13 voteeu
graph save g6, replace
reg mmx_progayindex2 c.eth_frac13 ethnocentric if sarajevocity==0 & yr2021==1, robust
mcp eth_frac13 ethnocentric
graph save g7, replace
reg mmx_progayindex2 c.eth_frac13 religious if sarajevocity==0 & yr2021==1, robust
mcp eth_frac13 religious
graph save g8, replace
graph combine g5.gph g6.gph g7.gph g8.gph
*Note: Additional Labeling and Formatting Required






*Online Appendix Material Replication

*Comparing Effects from the Sarajevo Pride Inside vs. Outside Sarajevo

*Figure 3:  Support for Sarajevo Pride In Sarajevo Over Time (%)
tab supportpride if sarajevocity==1 & treatment2==0
tab supportpride if sarajevocity==1 & treatment2==1

*Figure 4:  Support for Sarajevo Pride outside Sarajevo Over Time (%)
tab supportpride if sarajevocity==0 & treatment2==0
tab supportpride if sarajevocity==0 & treatment2==1
tab supportpride if sarajevocity==0 & treatment2==2

*Summary of Variables (pooled data, outside Sarajevo)
sum supportpride closetogay gayrights dgprogay progayindex2 i.treatment2 eth_frac13 diversity hadcovid diedcovid vaccinated ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0 & progayindex2~=.

*LGBT+ Support Index Components (Outside Sarajevo)

*Support for the Sarajevo Pride Over Time
tab supportpride if sarajevocity==0 & treatment2==0
tab supportpride if sarajevocity==0 & treatment2==1
tab supportpride if sarajevocity==0 & treatment2==2
*(Note: combine strongly and somewhat support/oppose responses)

*Social Distance to LGBT+ People Over Time
tab closetogay if sarajevocity==0 & treatment2==0
tab closetogay if sarajevocity==0 & treatment2==1
tab closetogay if sarajevocity==0 & treatment2==2
*(Note: record “not close at all” responses)

*Opposition to LGBT+ Rights Over Time
tab gayrights if sarajevocity==0 & treatment2==0
tab gayrights if sarajevocity==0 & treatment2==1
tab gayrights if sarajevocity==0 & treatment2==2
*(Note: record “strongly disagree” responses)

*Resource Allocation to Pro LGBT+ Group Over Time 
tab dgprogay if sarajevocity==0 & treatment2==0
tab dgprogay if sarajevocity==0 & treatment2==1
tab dgprogay if sarajevocity==0 & treatment2==2
*(Note: summarize responses <500, 500, and >500 MK)

*Figure 9:
twoway histogram supportpride if sarajevocity==0, percent discrete by(treatment2)

*Figure 10:
twoway histogram closetogay if sarajevocity==0, percent by(treatment2)

*Figure 11: 
twoway histogram gayrights if sarajevocity==0, percent by(treatment2)

*Figure 12:
twoway histogram dgprogay if sarajevocity==0, percent by(treatment2)

*Factor Analysis of LGBT+ Support Index (Dependent Variable)
factor supportpride closetogay gayrights dgprogay if sarajevocity==0
*use predict command to generate dgprogayindex2 (aldready generated)
*Note(we create a similar index called dgprogayindex which includes data from Sarajevo)
*use the 'norm' module package by Ansari and Mussida in STATA to normalize the index with the command 'norm progayindex2, method(mmx)' (already generated)

*LGBT+ Support Kernel Density Plot
kdensity progayindex2 if sarajevocity==0

*Manuscript Table 1 Robustness Checks

*Table 3:  Change in LGBT+ Support Over Time Outside Sarajevo (OLS Regression)
reg mmx_progayindex2 i.treatment2 if sarajevocity==0, cluster(municipality)
reg mmx_progayindex2 i.treatment2##c.eth_frac13 voteeu ethnocentric religious  if sarajevocity==0, cluster(municipality)
reg mmx_progayindex2 i.treatment2##c.eth_frac13  voteeu ethnocentric religious hadcovid diedcovid vaccinated Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0, cluster(municipality)

*Table 4. Change in LGBT+ Support Over Time Outside Sarajevo (Poisson Regression)
poisson mmx_progayindex2 i.treatment2 if sarajevocity==0, robust
poisson mmx_progayindex2 i.treatment2##c.eth_frac13 voteeu ethnocentric religious if sarajevocity==0, robust
poisson mmx_progayindex2 i.treatment2##c.eth_frac13 voteeu ethnocentric religious  hadcovid diedcovid vaccinated  Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0, robust

*Table 5. Change in LGBT+ Support Over Time Outside Sarajevo (Fixed Effects, OLS)
xtset municipality
xtreg progayindex2 i.treatment2  if sarajevocity==0, fe robust
xtreg progayindex2 i.treatment2  hadcovid diedcovid vaccinated ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0, fe robust

*Table 6. Change in LGBT+ Support Over Time Outside Sarajevo (Fixed Effects, Poisson Regression)
xtpoisson mmx_progayindex2 i.treatment2  if sarajevocity==0, fe robust
xtpoisson mmx_progayindex2 i.treatment2  hadcovid diedcovid vaccinated ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0, fe robust

*Table 7:  Change in LGBT+ Support Over Time Outside Sarajevo (OLS Regression)
reg progayindex2 i.treatment2 if sarajevocity==0, cluster(municipality)
reg progayindex2 i.treatment2##c.eth_frac13 voteeu ethnocentric religious  if sarajevocity==0, cluster(municipality)
reg progayindex2 i.treatment2##c.eth_frac13 voteeu ethnocentric religious hadcovid diedcovid vaccinated Bosniak Croat Serb Men age Education rural unemployed  rs log_casualty refugees ethnic_vote_share if sarajevocity==0, cluster(municipality)
*(note: response comparable with normalized dependent variable) 

*Table 8a:  Change in LGBT+ Support Over Time Outside Sarajevo (OLS Regression)
reg supportpride i.treatment2   if sarajevocity==0, robust
reg closetogay i.treatment2   if sarajevocity==0, robust
reg gayrights i.treatment2   if sarajevocity==0, robust
reg dgprogay i.treatment2   if sarajevocity==0, robust

*Table 8b:  Change in LGBT+ Support Over Time Outside Sarajevo (Ordered Probit Regression)
oprobit supportpride i.treatment2 if sarajevocity==0, robust
tobit closetogay i.treatment2 if sarajevocity==0, ll ul vce(robust)
tobit gayrights i.treatment2 if sarajevocity==0, ll ul vce(robust)
tobit gayrights i.treatment2 if sarajevocity==0, ll ul vce(robust)

*Table 9:  Covariable Balance Pre-Pride vs. Post-Pride
mlogit treatment2 ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0, cluster(municipality) baseoutcome(0)

*Table 10. Covariable Balance Immediate Post-Pride vs. Two-Year Follow-up
mlogit treatment2 ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0, cluster(municipality) baseoutcome(1)

*Table 11. Time Effects on LGBT+ Support (Coarsened Exact Matching)

*Immediately Pre/Post-Pride Outside Sarajevo

*downloading cem command below
ssc install cem

cem rural religious if treatment2~=2 & sarajevocity==0, treatment(treatment2)
reg progayindex2 treatment2  if sarajevocity==0, cluster(municipality)
reg progayindex2 treatment2  if sarajevocity==0 [pweight=cem_weights], cluster(municipality)

*Immediately Before 2019 Pride to Two Years After Pride
cem rural Bosniak if treatment2~=1 & sarajevocity==0, treatment(treatment2)
reg progayindex2 treatment2  if sarajevocity==0, cluster(municipality)
reg progayindex2 treatment2  if sarajevocity==0 [pweight=cem_weights], cluster(municipality)

*Immediately After 2019 Pride to Two Years After Pride
cem Bosniak if treatment2~=0 & sarajevocity==0, treatment(treatment2)
reg progayindex2 treatment2  if sarajevocity==0 & treatment~=0, cluster(municipality)
reg progayindex2 treatment2  if sarajevocity==0 & treatment~=0 [pweight=cem_weights], cluster(municipality)

*Table 12. Time Effects on LGBT+ Support (Propensity Score Matching)

*Immediately Before 2019 Pride to Two Years After Pride
teffects ipwra (progayindex2  ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs , linear) (treatment2 Bosniak Croat Serb Men age Education rural unemployed, probit) if treatment2~=1 & sarajevocity==0
teffects psmatch (progayindex2) (treatment2 ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs , probit) if treatment2~=1 & sarajevocity==0

*Immediately After 2019 Pride to Two Years After Pride
teffects ipwra (progayindex2  ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs , linear) (treatment2 Bosniak Croat Serb Men age Education rural unemployed, probit) if treatment2~=0 & sarajevocity==0
teffects psmatch (progayindex2) (treatment2 ethnocentric religious voteeu Bosniak Croat Serb Men age Education rural unemployed  rs , probit) if treatment2~=0 & sarajevocity==0

*Table 13:  Predictors of Covid-19 Loss, Infection, and Vaccination (Logit Regression)
logit diedcovid eth_frac13 religious ethnocentric voteeu hadcovid vaccinated Bosniak Croat Serb Men age Education rural unemployed rs if sarajevocity==0 & mmx_progayindex2!=. & treatment2==2, cluster(municipality)
logit hadcovid eth_frac13 religious ethnocentric voteeu diedcovid vaccinated Bosniak Croat Serb Men age Education rural unemployed rs if sarajevocity==0 & mmx_progayindex2!=. & treatment2==2, cluster(municipality)
logit vaccinated eth_frac13 religious ethnocentric voteeu hadcovid diedcovid Bosniak Croat Serb Men age Education rural unemployed rs if sarajevocity==0 & mmx_progayindex2!=. & treatment2==2, cluster(municipality)

*Table 14:  Cultural Receptivity Measures, Pearson’s Correlation Coefficient
pwcorr voteeu ethnocentric religious if sarajevocity==0, sig

*Table 15:  Cultural Receptivity Factor Analysis
factor eth_frac13 voteeu ethnocentric religious if sarajevocity==0

*Table 16:  Variance Inflation Factor
reg progayindex2 i.treatment2 c.eth_frac13  voteeu ethnocentric religious if sarajevocity==0, cluster(municipality)
vif

*Ethnic Fractionalization Histogram"
histogram eth_frac13 , normal percent

*Note: Code for alternative survey based "diversity" index provided in "Social Forces Replication Instructions". Variable already generated.

*Table 17:  Correlations with Individual Diversity Experience (OLS Regression)
reg diversity  voteeu ethnocentric religious if sarajevocity==0, cluster(municipality)

*Table 18:  Change in LGBT+ Support Over Time Outside Sarajevo (OLS Regression)
reg mmx_progayindex2 i.treatment2 diversity if sarajevocity==0, cluster(municipality)
reg mmx_progayindex2 i.treatment2##c.diversity if sarajevocity==0, cluster(municipality)
reg mmx_progayindex2 i.treatment2##c.diversity  voteeu ethnocentric religious hadcovid diedcovid vaccinated Bosniak Croat Serb Men age Education rural unemployed  rs if sarajevocity==0, cluster(municipality)

log close
